First report of natural infection of Anopheles gambiae s.s. and Anopheles coluzzii by Wolbachia and Microsporidia in Benin: a cross-sectional study.

BACKGROUND: Recently, bacterial endosymbiont, including Wolbachia and Microsporidia were found to limit the infection of Anopheles mosquitoes with Plasmodium falciparum. This study aimed to investigate the natural presence of key transmission-blocking endosymbionts in Anopheles gambiae and Anopheles coluzzii in Southern Benin. METHODS: The present study was conducted in seven communes (Cotonou, Porto-Novo, Aguégués, Ifangni, Pobè Athiémé, and Grand-Popo) of Southern Benin. Anopheles were collected using indoor/outdoor Human Landing Catches (HLCs) and Pyrethrum Spray Catches (PSCs). Following morphological identification, PCR was used to identify An. gambiae sensu lato (s.l.) to species level and to screen for the presence of both Wolbachia and Microsporidia. Plasmodium falciparum sporozoite infection was also assessed using ELISA. RESULTS: Overall, species composition in An. gambiae s.l. was 53.7% An. coluzzii, while the remainder was An. gambiae sensu stricto (s.s.). Combined data of the two sampling techniques revealed a mean infection prevalence with Wolbachia of 5.1% (95% CI 0.90-18.6) and 1.3% (95% CI 0.07-7.8) in An. gambiae s.s. and An. coluzzii, respectively. The mean infection prevalence with Microsporidia was 41.0% (95% CI 25.9-57.8) for An. gambiae s.s. and 57.0% (95% CI 45.4-67.9) for An. coluzzii. Wolbachia was only observed in Ifangni, Pobè, and Cotonou, while Microsporidia was detected in all study communes. Aggregated data for HLCs and PSCs showed a sporozoite rate (SR) of 0.80% (95% CI 0.09-2.87) and 0.69% (95% CI 0.09-2.87) for An. gambiae and An. coluzzii, respectively, with a mean of 0.74% (95% CI 0.20-1.90). Of the four individual mosquitoes which harboured P. falciparum, none were also infected with Wolbachia and one contained Microsporidia. CONCLUSIONS: The present study is the first report of natural infections of field-collected An. gambiae s.l. populations from Benin with Wolbachia and Microsporidia. Sustained efforts should be made to widen the spectrum of bacteria identified in mosquitoes, with the potential to develop endosymbiont-based control tools; such interventions could be the game-changer in the control of malaria and arboviral disease transmission.

Evidence of Transmission of Plasmodium vivax 210 and Plasmodium vivax 247 by Anopheles gambiae and An. coluzzii, Major Malaria Vectors in Benin/West Africa.

Current diagnostic and surveillance systems in Benin are not designed to accurately identify or report non-Plasmodium falciparum (Pf) human malaria infections. This study aims to assess and compare the prevalence of circumsporozoite protein (CSP) antibodies of Pf and P. vivax (Pv) in Anopheles gambiae s.l. in Benin. For that, mosquito collections were performed through human landing catches (HLC) and pyrethrum spray catches (PSC). The collected mosquitoes were morphologically identified, and Pf, Pv 210, and Pv 247 CSP antibodies were sought in An. gambiae s.l. through the ELISA and polymerase chain reaction (PCR) techniques. Of the 32,773 collected mosquitoes, 20.9% were An. gambiae s.l., 3.9% An. funestus gr., and 0.6% An. nili gr. In An. gambiae s.l., the sporozoite rate was 2.6% (95% CI: 2.1-3.1) for Pf, against 0.30% (95% CI: 0.1-0.5) and 0.2% (95% CI: 0.1-0.4), respectively, for Pv 210 and Pv 247. P. falciparum sporozoite positive mosquitoes were mostly An. gambiae (64.35%), followed by An. coluzzii (34.78%) and An. arabiensis (0.86%). At the opposite, for the Pv 210 sporozoite-positive mosquitoes, An. coluzzii and An. gambiae accounted for 76.92% and 23.08%, respectively. Overall, the present study shows that P. falciparum is not the only Plasmodium species involved in malaria cases in Benin.

Pesticides and the evolution of the genetic structure of Anopheles coluzzii populations in some localities in Benin (West Africa).

BACKGROUND: Changes in the natural habitats of insect groups are determined the genetic polymorphisms between individuals. The objective of this study was to establish the genetic structure of the Anopheles coluzzii populations in four localities of Benin. METHODS: Insecticide surveys and larval sampling were conducted on 4 study localities, including Cotonou, Ketou, Zagnanado, and Sô-Ava. Molecular characterizations were performed on the Anopheles mosquitoes collected with the allelic and genotypic frequencies of kdr gene determined. The multiple comparison Chi square test for proportions was performed with R version 3.3.3. Next, the observed heterozygosity, expected heterozygosity, and indices of fixation, and genetic differentiation were estimated. Finally, the Hardy-Weinberg equilibrium (EHW) was determined to assess whether panmixia exists in the different populations of mosquitoes of the agroecological zones under study. RESULTS: Carbamates, pyrethroids, organophosphorus and organochlorines use have been reported in all localities except Sô-Ava. Anopheles coluzzii was strongly represented across all study localities. The L1014F allele was observed in the localities of Kétou, Cotonou and Zagnanado. Likewise, insecticide selection pressure of homozygous resistant individuals (L1014F/L1014F) was significantly higher in Kétou, Cotonou and Zagnanado (p value < 0.05). Surprisingly in Sô-Ava, a relatively high frequency of the L1014F allele despite the reported absence of pesticide use was observed. All mosquito populations were found to be deficient in heterozygosity across the study sites (FIS< 0). No genetic differentiation (FST< 0) was observed in the localities of Zagnanado and Kétou. CONCLUSION: The survey on the use of insecticides showed that insecticide selection pressures differ across the investigated localities. It would be desirable to rotate or apply formulations of combined products with different modes of action. Doing so would enable a better management of resistant homozygous individuals, and mitigate the resistance effect of commonly used insecticides.

The current distribution and characterization of the L1014F resistance allele of the kdr gene in three malaria vectors (Anopheles gambiae, Anopheles coluzzii, Anopheles arabiensis) in Benin (West Africa).

BACKGROUND: The fight against malaria faces various biological obstacles, including the resistance of parasites to anti-malarial drugs and the resistance of mosquito vectors to insecticides. The resistance of Anopheles gambiae sensu lato (s.l.) to pyrethroids, the only class of insecticides used to impregnate mosquito nets, is known in Benin; the expansion of this resistance is influenced by the existence of gene flow between species, otherwise by the presence or absence of the kdr mutation in them. The objective of this study is to determine the spatial distribution of An. gambiae and the level of expression of the pyrethroid resistance kdr gene in seven agro-ecological zones of Benin. METHODS: The study was conducted in 18 localities belonging to seven agro-ecological zones where environmental parameters varied. The sites represent the main areas of eco-epidemiological malaria in Benin. Anopheles gambiae larvae were collected in natural breeding sites using ladles and dipping method and reared under standard conditions. These larvae were reared under standard conditions of temperature and humidity (26 to 30 °C and 60 to 90%) at the insectarium of the Centre de Recherche Entomologique de Cotonou (CREC). Adult female mosquitoes having emerged are morphologically and molecularly identified. Homozygous resistant (1014F/1014F), homozygous sensitive (1014L/1014L) and heterozygous (1014F/1014L) genotypes of the L1014F kdr gene mutation are determined by PCR. RESULTS: A total of 677 An. gambiae was subjected at the PCR. The results revealed the presence of three vector species of the An. gambiae complex, of which 409 Anopheles coluzzii, 259 An. gambiae, 5 hybrids (An. coluzzii/An. gambiae) and 4 Anopheles arabiensis in the different agro-ecological zones. The four An. arabiensis were only found in Dassa, a locality in the cotton zone of central Benin. The frequency of distribution of the L1014F allele of the kdr gene varies from 84.48 to 100% in An. gambiae, from 80 to 100% in An. coluzzii and from 0 to 75% in An. arabiensis in the different agro-ecological zones. Moreover, a significant difference is generally observed in the distribution of the L1014F allele (P < 0.05). By comparing in pairs the distribution frequencies of this allele in the two species by agro-ecological zone, only a significant difference is noted in the central cotton and fishery zones (P = 0.0496). CONCLUSION: In summary, even if the data are in small portions, the An. Arabiensis species was found only in central Benin and the L1014F allele of the kdr gene is widespread and seems to fix in all the species recorded in the different agro-ecological zones. This situation amplifies the problem of resistance, which could eventually be a significant obstacle for the malaria vectors control. Similarly, a study of their genetic structure via the L1014F allele is necessary in order to put in place strategies to manage this resistance. These strategies will take into account both the ecology and the genetic diversity of the organisms involved to preserve the effectiveness of pyrethroids, the only insecticides used for the impregnation of mosquito nets.

Involvement of Anopheles nili in Plasmodium falciparum transmission in North Benin.

BACKGROUND: Several studies carried out in Benin have shown the involvement of Anopheles gambiae sensu lato (s.l.), the Anopheles funestus group in malaria transmission, but none of them reported the contribution of the Anopheles nili group to the transmission of this disease. The current study investigated the question through an entomological cross-sectional survey performed in Northern Benin. METHODS: Mosquito samplings were performed in September and October 2017 in 4 villages located in two districts: Bambaba and Wodara (Kérou district) and, Péhunco 2 and Béké (Péhunco district). The collections were carried out indoors and outdoors using human landing catches (HLC) to assess the human biting rate (HBR) and pyrethrum spray catches (PSC) to evaluate the blood feeding rate and the blood meal origin using the ELISA test. All collected mosquitoes were morphologically identified and, the polymerase chain reaction (PCR) technique was used for molecular identification of sibling species of An. gambiae s.l., An. funestus group and An. nili group sporozoite index (SI) was also assessed by the ELISA test. RESULTS: Overall, An. gambiae s.l., An. funestus group and An. nili group were the three vectors found in the study area. A significantly higher human biting rate (HBR) was recorded in An. nili group (5 bites/human/night) compared to An. funestus group (0.656 bites/human/night) in the Kérou district (p < 0.0001). Anopheles gambiae s.l. displayed the highest HBR (26.19 bites/human/night) in the same district. The entomological inoculation rate (EIR) was 1.875 infected bites/human/month in An. nili group against 13.05 infected bites/human/month in An. gambiae s.l. and 0.938 infected bites/human/month in An. funestus group in Kérou. In Péhunco, the EIR was 1.02 infected bites/human/month in An. gambiae s.l. PCR results showed that An. nili sensu stricto (s.s.) and An. funestus s.s. were the only species of the An. nili and An. funestus groups, respectively. The anthropophagic character of An. gambiae s.l. was also highlighted. CONCLUSION: This study provides useful information on the contribution of An. nili group as secondary vector to malaria transmission in northern Benin. Broader studies must also be carried out in a larger study area to assess the involvement of other Anopheles species to malaria transmission. This will aid to better plan malaria vector control interventions.

Genetic structure of Anopheles gambiae s.s populations following the use of insecticides on several consecutive years in southern Benin.

BACKGROUND: Several studies have reported the strong resistance of Anopheles gambiae s.l. complex species to pyrethroids. The voltage-dependent sodium channel (Vgsc) gene is the main target of pyrethroids and DDT. In Benin, the frequency of the resistant allele (L1014F) of this gene varies along the north-south transect. Monitoring the evolution of resistance is necessary to better appreciate the genetic structure of vector populations in localities subject to the intensive use of chemicals associated with other control initiatives. The purpose of this study was to map the distribution of pyrethroid insecticide resistance alleles of the Kdr gene in malaria vectors in different regions and ecological facies in order to identify the evolutionary forces that might be the basis of anopheline population dynamics. METHODS: The characterization of Anopheles gambiae s.l. populations and resistance mechanisms were performed using adult mosquitoes obtained from larvae collected in the four agroecological zones in southern Benin. Genomic DNA extraction was performed on whole mosquitoes.The extracted genomic DNA from them were used for the molecular identification of species in Anopheles gambiae s.l. complex and the identification of genotypes related to pyrethroid resistance as the Kdr gene amino acid position 1014 in sodium channel. Molecular speciation and genotyping of Kdr resistant alleles (1014) were done using PCR.Genepop software version 4.2 was used to calculate allelic and genotypic frequencies in each agroecological zone. The p value of the allelic frequency was determined using the binomial test function in R version 3.3.3. The Hardy-Weinberg equilibrium was checked for each population with Genetics software version 1.3.8.1. The observed heterozygosity and the expected heterozygosity as well as the fixation index and genetic differentiation index within and between populations were calculated using Genepop software version 4.2. RESULTS: During the study period, Anopheles coluzzii was the major species in all agroecological zones while Anopheles gambiae was scarcely represented. Regardless of the species, resistant homozygote individuals (L1014F/L1014F) were dominant in all agroecological zones, showing a strong selection of the resistant allele (L1014F). All populations showed a deficit of heterozygosity. No genetic differentiation was observed between the different populations of the two species. For Anopheles coluzzii, there was a small differentiation among the populations of the central cotton and bar-lands zones. The genetic differentiation was modest among the population of the fisheries zone (Fst = 0.1295). The genetic differentiation was very high in the population of Anopheles gambiae of the bar-lands zone (Fst = 0.2408). CONCLUSION: This study revealed that the use of insecticides in Benin for years has altered the genetic structure of Anopheles gambiae s.s. populations in all agroecological zones of southern Benin. It would be desirable to orientate vector control efforts towards the use of insecticides other than pyrethroids and DDT or combinations of insecticides with different modes of action.

Evidence of multiple insecticide resistance mechanisms in Anopheles gambiae populations in Bangui, Central African Republic.

BACKGROUND: Knowledge of insecticide resistance status in the main malaria vectors is an essential component of effective malaria vector control. This study presents the first evaluation of the status of insecticide resistance in Anopheles gambiae populations from Bangui, the Central African Republic. METHODS: Anopheles mosquitoes were reared from larvae collected in seven districts of Bangui between September to November 2014. The World Health Organisation's bioassay susceptibility tests to lambda-cyhalothrin (0.05%), deltamethrin (0.05%), DDT (4%), malathion (5%), fenitrothion (1%) and bendiocarb (0.1%) were performed on adult females. Species and molecular forms as well as the presence of L1014F kdr and Ace-1 R mutations were assessed by PCR. Additional tests were conducted to assess metabolic resistance status. RESULTS: After 1 h exposure, a significant difference of knockdown effect was observed between districts in all insecticides tested except deltamethrin and malathion. The mortality rate (MR) of pyrethroids group ranging from 27% (CI: 19-37.5) in Petevo to 86% (CI: 77.6-92.1) in Gbanikola; while for DDT, MR ranged from 5% (CI: 1.6-11.3) in Centre-ville to 39% (CI: 29.4-49.3) in Ouango. For the organophosphate group a MR of 100% was observed in all districts except Gbanikola where a MR of 96% (CI: 90-98.9) was recorded. The mortality induced by bendiocarb was very heterogeneous, ranging from 75% (CI: 62.8-82.8) in Yapele to 99% (CI: 84.5-100) in Centre-ville. A high level of kdr-w (L1014F) frequency was observed in all districts ranging from 93 to 100%; however, no kdr-e (L1014S) and Ace-1 R mutation were found in all tested mosquitoes. Data of biochemical analysis showed significant overexpression activities of cytochrome P450, GST and esterases in Gbanikola and Yapele (χ 2 = 31.85, df = 2, P < 0.001). By contrast, esterases activities using α and ß-naphthyl acetate were significantly low in mosquitoes from PK10 and Ouango in comparison to Kisumu strain (χ 2 = 17.34, df = 2, P < 0.005). CONCLUSIONS: Evidence of resistance to DDT and pyrethroids as well as precocious emergence of resistance to carbamates were detected among A. gambiae mosquitoes from Bangui, including target-site mutations and metabolic mechanisms. The co-existence of these resistance mechanisms in A. gambiae may be a serious obstacle for the future success of malaria control programmes in this region.

Six years of experience in entomological surveillance of indoor residual spraying against malaria transmission in Benin: lessons learned, challenges and outlooks.

BACKGROUND: From 2008 to 2013, a prevention intervention against malaria based on indoor residual spraying (IRS) was implemented in Benin. From 2008 to 2012, Ficam M(®), a bendiocarb-containing product was used for house spraying, in association with pirimiphos methyl EC (Actellic EC) in 2013. This operation aimed to strengthen the effectiveness of treated nets so as to expedite the achievement of Millennium Development Goals (MDGs): the reduction of morbidity and mortality due to malaria by 75 % from 2000 to 2015. METHODS: Monitoring and evaluation (M&E) was implemented in order to evaluate the impact of IRS intervention on malaria transmission. Anopheles gambiae s.l. populations were sampled by human landing catch. In addition, window exit traps and pyrethrum spray catches were performed to assess exophagic behaviour of Anopheles vectors the main malaria vector in the treated areas. The residual activity of insecticide in the treated walls was also assessed using WHO bioassay test. RESULTS: The purpose of this project was to draw attention to new challenges and future prospects for the success of IRS in Benin. The main strength of the intervention was a large-scale operation in which more than 80 % of the houses were treated due to the strong adhesion of population. In addition, a significant reduction of the EIR in areas under IRS were observed. However, there were many challenges including the high cost of IRS implementation and the identification of suitable areas to implement IRS. This was because of the low and short residual effect of the insecticides recommended for IRS and the management strategy for vector resistance to insecticides. This indicated that challenges are accompanied by suggested solutions. For the cost of IRS to be accessible to states, then local organizations need to be created in partnership with the National Malaria Control Programme (NMCP) in order to ensure relevant planning and implementation of IRS. CONCLUSION: As an anticipatory measure against vector resistance, this paper proposes various methods, such as periodic IRS based on a combination of two or three insecticides of different classes used in rotation every two or three years.

A shift from indoor residual spraying (IRS) with bendiocarb to long-lasting insecticidal (mosquito) nets (LLINs) associated with changes in malaria transmission indicators in pyrethroid resistance areas in Benin.

BACKGROUND: Indoor residual spraying (IRS) was implemented in the department of Ouémé-Plateau, southern Benin, in 2008 and withdrawn in 2011, when long lasting insecticidal nets (LLINs) were distributed to the communities that were previously targeted by IRS. Did the LLIN strategy provide a better level of protection against malaria transmission than IRS? METHODS: Entomological surveillance was carried out to assess indicators of transmission risk during the last year of IRS and the first year after the LLIN intervention was put in place (2010-2011). Mosquito biting density was sampled by human landing collection (HLC). Females of Anopheles gambiae s.l. were dissected to estimate the parity rates and the blood meal index. A subsample of the An. gambiae s.l. collection was tested for presence of Plasmodium falciparum sporozoites. In addition, window exit traps and pyrethrum spray catches were performed to assess exophagic behavior of Anopheles vectors. RESULTS: There were significant increases in all the indicators following withdrawal of IRS. Vector biting density (p<0.001) and longevity (OR=3.81[3.01-4.82] 95% CI; p<0.001) of the An. gambiae s.l. increased significantly; so too did the blood meal index (OR=1.48 [1.1-1.99] 95% CI; p<0.001). Entomological inoculation rate, after IRS withdrawal at one surveillance site, Adjohoun, rose two fold (9.0 infected bites/person/9 months (Apr-Dec 2011) versus 3.66 infective bites/person during the 9 months preceding IRS (Apr-Dec 2010). A second site, Missérété, experienced a six-fold increase after IRS cessation (15.1 infective bites/person/9 months versus 2.41 during IRS). Exophily after IRS cessation decreased significantly in all areas (p<0.001) suggesting that mosquitoes were more likely to rest in houses with LLINs, than in houses subjected to IRS. CONCLUSION: LLINs did not impact on indicators of transmission to the same levels as did IRS after IRS withdrawal.